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Observation of tree-level B decays with ss Production from Gluon Radiation 
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We report on our search for ^ oi*''^ K' n' , B° Di*>~^ K'^tt' , and ^ Dr'^K-R- de- 
cays in 383 million T(45') ^ -B-B events collected by the BABAR detector at the PEP-II asymmetric- 
energy B-factory. The decays proceed via a tree-level b —f c quark transition in which a gluon 
radiates into an ss pair. Their branching fractions are measured to be B{B~ Df K~ n~) = 
(2.02 ± O.Ustat ± O.SSsyst) X 10"*, B{B- D*+K~-k-) = (1.67 ± Q.lQstat ± 0.35,j,,t) x 10"*, 
B{B° -» D+KItt-) = (0.55±0.13,tat±0.10,«s^)xlO-^ and B(B- -» D+ K') = (0.11±0.0W± 
0.02syst) x 10"". Upper limits at the 90% C.L. are set on B{B° -* D'+Kli^-) < 0.55 x 10"* and 
B(B~' Dl'^K~K~) < 0.15 X 10~*. We present evidence that the invariant mass distributions 
of _Di*'^_ft'~ pairs from B~ — > D'f''^K~n~ decays are inconsistent with the phase-space model, 
suggesting the presence of charm resonances lying below the oi*^'^ K~ threshold. 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 

Evidence for inclusive flavor correlated production of 
Z?^ in B~ decays was reported recently [l| with a branch- 
ing fraction of 6(5" D+X) = (1.2 ±0.4)% 0]. Along 
with D*^X decays, these decays are mediated by 

a b c quark transition and require at least three fi- 
nal state particles, including the production of an ss pair 
from the vacuum via radiative gluon pair production. Ex- 
amples for three-body B~ decays with a Ds in the fi- 
nal state are B~ Ds K~'k~ . The Feynman diagram 
for B- Di*^^K'n- decays are shown in Figure [1] 
The corresponding decays are B'^ I?s*''^A'°7r". 
By replacing the tt~ in Figure [T] with a , we get the 
Cabibbo-suppressed decays B~ D^*^^ K~ K~ . 




FIG. 1: Feynman diagram for B Dl K n . 

In addition to the dominant diagram, B^ 
oi*''^ K^TT^ can occur via the color-suppressed diagram 
where the constituent u's of the K~ and tt~ are switched. 
Although a color-suppressed contribution does not exist 
for B° Di*'^^K° TT , a sub-dominant contribution from 
a VF-exchange diagram with ss and dd popping may ex- 
ist instead. Either of these contributions could cause a 
deviation from the naive expectation of two for the ratio 
of B- Di*^+K-TT- to B" ^ Di*^^K°Tr- branching 
fractions. 

The Ds -K pair could come from intermediate charm 
resonances instead of directly from the B. It has been 
proposed that these resonances can play a significant role 
in B^ DfK~Tr~ decays 0] despite their masses lying 
below the m{DsK) production threshold [3l- In this case, 
it may be possible to measure the parameters of the res- 



onances such as their masses and widths, complementary 
to the analysis using B Dnn decays 

Along with exclusive B^ — > Z?i*-'^X and 
Ds X three-body decays, no decays proceeding via ra- 
diative gluon ss pair production at the tree level have 
hitherto been observed. Upper limits on the branch- 
ing fractions of the B" —> Di*''^ K^tt' and B" 
Di*^^K^n- modes have been placed by ARGUS Q. In 
this paper we report first observations of the decay modes 
B~ Di*^^K-n- and B^ D+K°t:-, evidence for 
B^ Df K~ , and limits on the branching fractions 
of ^0 -> D*+K°n- and B" ^ Dl+R-R-. We also 
present a measurement of an invariant mass distribution 
of D^s^^ K~ pairs from B^ Di*'^/! "tt" decays and 
compare it to the spectrum obtained from a phase-space 
model. 

The analysis uses approximately 383 million TiAS) — > 
BB events created by the PEP-II e+e" collider and col- 
lected by the BaBAR detector. The BAR4R. detector is 
described elsewhere @. 

Optimal selection criteria and probability density func- 
tions of selection variables are determined by an analysis 
based on Monte Carlo (MC) simulation of both signal 
and background events. We use GEANT4 Q software 
to simulate interactions of particles traversing the BaBAR 
detector, taking into account the varying detector condi- 
tions and beam backgrounds. We verify with MC sim- 
ulation that resolutions and background levels correctly 
describe the data. 

Candidates for 13+ mesons are reconstructed in the 
modes D+_ ^ (j)TT+ , K*°K+, and K°K+, with ^ 
K+K-, K*° -> K-TT+ and K^^ -* tt+tt". The can- 
didates are reconstructed from two oppositely-charged 
tracks that come from a common vertex displaced from 
the e^e~ interaction point. We require the significance 
of this displacement (the measured iCg flight distance 
divided by its estimated error) to exceed 2. All other 
tracks are required to originate less than 1.5 cm away 
from the e+e~ interaction point in the transverse plane 
and less than 10 cm along the beam axis. Charged kaon 
candidates must satisfy identification criteria that are 
typically around 92% efficient [§| , depending on momen- 
tum and polar angle, and have a pion misidentification 
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TABLE L Summary of results for the total detection efficiencies e excluding the subsequent branching fractions of Di*^ decay 
modes {D*^ Df'y, — > (?!)7r+, K'^K'^ , K^K'^), expected peaking background events ripeaking with statistical uncertainties 
from fits of the ities distributions obtained using the invariant mass sidebands, final signal (risig) and background (ribkg) 
yields with statistical uncertainties from niES fits adjusted to account for estimated peaking backgrounds, and cross-feed 
contributions, branching fractions B with statistical and systematic uncertainties, significances s{a) calculated by comparing 
the likelihood maximum of the nominal fit to that of the fit with the signal yield fixed to the difi'erence between the raw and 
corrected signal yields, and upper limits UL on the branching fractions for B° D*+K'^n- and ^ D'+K-R-. The 
background yields ntkg are calculated in the region 5.27 — 5.29 GeV/c^ . 

Mode e^^ £k*k Sk°k 'T-pcaking Usig ribkg B X 10"" s((j) UL (90% C.L.) 

B-^DtK-jT- 11.1% 6.8% 9.6% 41 ± 9 430 ± 29 182 ± 6 2.02 ± 0.13 ± 0.38 21 — 
B-^D*+K-TT- 5.9% 3.6% 5.1% 4±5 178 ± 18 87.1 ±3.5 1.67 ± 0.16 ± 0.35 14 — 
B^^D^ifsTT" 8.8% 5.3% 7.6% 28 ± 6 61.8 ± 14.4 94.5 ± 5.5 0.55 ± 0.13 ± 0.10 5.2 — 
B"^D;+K°^7v- 3.8% 2.3% 3.4% -1.1 ±2.7 13.6 ± 8.4 62.8 ± 3.4 0.29 ± 0.18 ± 0.073 1.8 0.55 x 10"* 
B-^DtK-R- 7.1% 4.3% 6.3% -0.3 ±1.9 14.4 ± 5.6 9.8 ± 1.3 0.11 ± 0.04 ± 0.02 3.3 — 
B- Dt"^ R- R- 3.8% 2.4% 3.5% -1.7 ±1.3 4.7 ±4.0 6.5 ± 0.9 0.07 ± 0.06 ± 0.02 1.3 



rate at the 5% level. The cj) K+K', K*° K--K+ 
and — > tt+tt" candidates are required to have in- 
variant masses within ±15 McV/c^, ±50 MeV/c^ and 
±10 MeV/c^ of their nominal masses, respectively 
The full polarization of the and (j) mesons from the 
D'^ decays are exploited to reject backgrounds through 
the use of the helicity angle Oh, defined as the angle 
between the momentum vector and the direction of 
flight of the £)+ in the K*'^ or cj) rest frame. The K*^ and 
(j) candidates are required to have | cos Oh \ greater than 
0.5. 

The candidates arc reconstructed in the mode 

Z?*+ Dfj. The photons are accepted if their energy is 
greater than 100 MeV. Photons from D*~^ candidates are 
rejected if, when combined with any other photon having 
an energy greater than 150 MeV, they belong to a pho- 
ton pair whose invariant mass lies within ±10 MeV/c^ of 
the 7r° mass. The Df candidates are required to have 
invariant masses in the interval ±10 MeV/c^ of the nom- 
inal £)+ mass while the invariant masses of can- 
didates he in the range from m(D*+) — 15 MeV/c^ to 
m{D*+) + 10 McV/c^. All D+ candidates are subjected 
to a mass-constrained fit after selection. The invariant 
mass of the D*"'" is calculated after the mass constraint 
on the daughter D+ has been applied. Subsequently, all 
D*~^ candidates are subjected to mass-constrained fits. 
To eliminate B° oi*^^ D- , D- K°tt- events from 
the 5° Di*^~^ K'^tt'^ samples, the invariant mass of the 
K'^ and 7r~ must be outside a 40 MeV/c^ window around 
the D~ mass. 

Finally, the B meson candidates are formed using the 
reconstructed combinations of D^K~tt~ , _D*+A'~7r~, 
D+K°Tr-, D*+K°Tr-, D+K~K', and D*+K'K~. 

The background from continuum qq production (where 



q = u,d,s,c) is suppressed based on the event topology. 
The event shape variables, i?2 (the ratio of the second to 
zeroth Fox- Wolfram moments [10|) and L2/L0 (the ratio 
of the second and zeroth angular moments of the energy 
flow about the B thrust axis [ll[), are combined in a 
Fisher discriminant (J-) to effectively exploit the differ- 
ence between the shapes of e+e~ — s- BB and e+e" qq 
events. A selection is applied to J- such that 80% of 
continuum background is rejected while maintaining 80% 
signal efficiency. 

The signals are extracted using the energy-substituted 
mass mEs = \J — (^ - p,*)^ and the energy difference 
Ai? = Yl,% + P*^ - -E^b 1 where E'jJ is the beam en- 
ergy in the laboratory frame, p* is the momentum of the 
daughter particle i of the B meson candidate also in the 
laboratory frame, and mi is the mass hypothesis for par- 
ticle i. For signal events, toes peaks at the B meson mass 
with a resolution of about 2.6 MeV/ (? and Ai? peaks near 
zero with a resolution of 13 MeV. The B candidates are 
required to have |A£:| < 25 MeV and toes > 5.2 GeV/c^. 
After all selection criteria are applied, we find the frac- 
tion of events containing more than one B candidate to 
be between 3% and 11% depending on the decay mode. 
In these instances, the B candidate with IS.E closest to 
zero is chosen. The estimated B reconstruction efficien- 
cies, excluding the subsequent branching fractions, are 
shown in Table HI 

Background events that pass these selection criteria 
are represented by approximately equal amounts of qq 
continuum and BB events. We parametrize their toes 
distributions by a threshold function [l^ : 

/(toes) ^ TOES \/l-a?^exp[-^(l - 2^^)], 

where x = 2m-E.s/ y/s, y/s is the total energy of the beams 
in their center of mass frame, and ^ is a fit parameter. 



6 



A study using simulated B decays reveals significant 
numbers of background events peaking in the regions of 
5.272 < Toes < 5.288 GeV/c^ and |A£;| < 25 MeV sim- 
ilar to the reconstructed signal candidates. This peak- 
ing background is due to charmless and charmonium B 
decays with the same set of particles as signal in the 
final state. The peaking contribution is evaluated using 
the data by reconstructing Z?s ' K~ii~,Dl ' A^tt" and 
Di*^~^ K~ combinations, where "D^" candidates are 
selected from 25 - 40 MeV/c^ sidebands around the Df 
nominal mass. In this procedure, we use the same selec- 
tion requirements as for the signal except that "1)^" can- 
didates are not mass constrained. Studies revealed that 
constraining the Df mass did not significantly affect the 
resolutions of toes and AE distributions and that events 
in the Df mass sidebands are a good representation of 
the background under the Z)+ peak. Table |T] shows the 
fit yields of the peaking background contribution under 
the Toes peak for each mode. 

A matrix is constructed to study the cross-feed be- 
tween the signal modes. Its elements describe the contri- 
butions of each mode according to the levels seen in MC 
samples. No off-diagonal element of the cross-feed ma- 
trix exceeds 2%; this near-diagonal structure indicates 
effective suppression of the cross-feed contributions by 
application of the selection criteria. 

Figure [2] shows the toes spectra of the reconstructed 
B candidates. For each mode, we perform an extended 
unbinncd maximum likelihood (ML) fit to the toes distri- 
butions using the candidates from all Df decay modes 
combined. The toes distributions are fit with the sum 
of two functions: /(toes) characterizing the combinato- 
rial background and a Gaussian function to describe the 
signal. The likelihood function is given by: 
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FIG. 2: niES spectra for the B — > DfK tv (top left), 
B- D*+K-'K- (top right), B° -> D+A'^tt" (middle left), 
B° -» D*+Kl-K- (middle right), B" D+ R- K~ (bottom 
left), and B^ D^'^ K~ K~ (bottom right). Solid curves 
show the fit results, as explained in the text. Dashed lines in 
the signal regions correspond to the peaking and non-peaking 
background components of the fit. The data are the points 
with error bars. 



p-(ns>g+nbfcg) IL 
i=l 

where Pp^ and P/''^^ are the probability density func- 
tions for the signal and background, Ugig and ribkg are 
the number of signal and background events, and N is 
the total number of events in the fit 

The final signal yields are obtained by subtracting the 
estimated peaking background and cross-feed contribu- 
tions from the yields of the niES fits described in the pre- 
ceding paragraph. No peaking background is subtracted 
from modes that have ripcaking less than zero in Table |I] 
since these values are consistent with zero although their 
errors are still propagated. The final values are given in 
the Usig column of Table HI The total signal yield in each 
B decay mode is related to the B branching fraction B 
using the following expression: 

B = nsig/{NBB - ^B^-ei), 



where Ngj^ is the number of produced BB pairs, Bi is 
the product of the intermediate branching ratios, is the 
reconstruction efficiency (from Table |T| and the sum is 
over D+ modes (i = 07r+, K*°K+ , KgK+). As an input 
to the calculations, we used branching fraction numbers 
from 0. The results of these calculations are summa- 
rized in Table m 

The total relative systematic uncertainty in the B 
branching fractions is estimated to be approximately 
19% — 25% depending on the decay mode. The largest 
contribution, an uncertainty of 15%, comes from the Df 
branching fractions. The differences between selection 
efficiencies in MC and in the data (estimated using the 
control mode B^ D° , D° K-t:+) contribute to 

the systematic uncertainty (5% — 10%) as does the effi- 
ciency dependence on the Ds AT invariant mass spec- 
trum (7% — 9%). In the toes fits of the lower statistics 
modes (£)*+A'°7r~ and D*'^K~K~) the signal Gaussian 
parameters and ^/s in /(toes) are fixed to ensure fit con- 
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vergence. The associated systematic uncertainties are 
14% and 9%, respectively. The entries in the cross-feed 
matrix affecting the Di*^^ modes vary by 8%(5%) 
when they are calculated with MC events weighted ac- 
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cording to the observed spectra of the Dl 
mass. 

For the 5° ^ D*+/v!?7r- and B- D*+ R- R- de- 
cay modes, the upper limits arc set using a frequcntist 
approach [9| and taking into account the systematic un- 
certainties. The upper limits arc summarized in Table [H 

(*)+ — 

Studies of the invariant mass spectra of the Dl 
system in Dl K~'k~ modes reveal distributions 

incompatible with those of three-body phase space. As 
shown in Figure[3l there are enhancements in the number 
of events at the lower ends of the m{D^^'^~^ K~) spectra. 
These enhancements suggest the presence of charm reso- 
nances lying below the d'-*^^K- threshold 0. 

In summary, B~ DfK~iT~, B~ -> DI'^K'tt^ and 
DfKgTT^ decays are observed for the first time 
each with significance greater than Scr. Evidence for 
B~ ^ D^K~K~ was found with a significance slightly 
greater than 3(T. Upper limits are set on the branching 
fractions of the two decay modes with significances lower 
than 2a: B° D*+K°tt- and ^ Dl+K'K-. 

The ratios of B{B- D^*^^ R- R-) to B{B- 

TT ) are consistent with the expected Cabibbo 
suppression. The branching fraction of B^ DfK'^7T~ 
is less than half that of B^ K~tt~ ; this may be 

due to the W-exchangc diagram correction to the neu- 
tral mode and the color-suppressed contribution to the 
charged mode. 
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FIG. 3: Di*^+K- invariant mass spectra for the B — > 
DfK~n~ (left) and B^ — > D*^ tt^ (right) decay modes 
using the data. A requirement of m^;^ > 5.270 GeV/c^ is ap- 
plied to the events shown in the figure, in addition to the sig- 
nal selection described in the text. Combinatoric background 
is approximated and then subtracted using events outside the 
mBS signal region {tues < 5.265 GeV/c^). The histogram 
shows the non-resonant signal MC events distribution, scaled 
to the number of events in the data signal region. 



